Searches for CP, T and CPT violation in B°B° mixing at BABAR 
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We describe measurements of CP and CPT violation parameters in B°B° oscillations, performed 
at BABAR by using BB events collected at the T(4S) mass peak. These are obtained from two 
different ^-reconstruction methods, namely: searching for two inclusive semileptonic B° decays in 
an event or partially reconstructing one of the B° mesons in the D*~£ + vi channel. The results are 
in agreement with one another and with the most recent theoretical calculations in the Standard 
Model scenario. 
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I. THEORETICAL OVERVIEW 

CP violation due to the neutral- B meson oscillation 
processes can be parameterized in two different scenar- 
ios: either in the Standard Model (SM) field theory, with 
some extensions allowed within the CPT-conservation 
hypothesis, or, more generally, also taking into account 
the possibility of CPT violation. 

In the SM case, B°B° mixing is accounted for using 
a time-dependent perturbation theory that leads to an 
effective hamiltonian H e // = M — ~T, that is the 2x2 
projection of the actual interaction operator on the B° 
and B° states only. In solving the eigenvalue problem, 
linear relations are found connecting the B flavor eigen- 
states to the physical states: 
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In the Standard Model, the magnitude of the ratio q/p 
is very nearly unity: 
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not just because jP^I is small, but because the CP- 
violating quantity Im(Ti2/Mi2) is suppressed by an addi- 
tional factor (m^ — )/m 2 b w 0.1 relative to \Ti 2 /M 12 \. 
When the remaining factors are included, the expectation 
is |Im(r 12 /Af 12 )| < lO" 3 0,0. 

The parameter \q/p\ is sensitive to New Physics (NP) 
within a CPT-conserving theory Assuming 3- 

generation unitarity of the CKM matrix Q and tree-level 
processes dominated by the SM, the CP asymmetry can 
be modified as follows: 
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Cd and <fid being a generic NP amplitude and phase. 
Within the NP parameter space, \q/p\ — 1 could be en- 
hanced up to an order of magnitude. 



If CPT violation is also allowed, its effects are usually 
parameterized in terms of an additional complex param- 
eter z in the eigenstate definitions JTJ): 
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Kostelecky [5] first pointed out that such parameteriza- 
tion using a constant value of z is not completely satisfac- 
tory, since CPT violation in a field theory implies that 
Lorentz covariance, too, is violated. A more appropri- 
ate parameterization is achieved by setting z oc fl^Aa^, 
where (3^ = 7b(1,/3b) is the decaying f?-meson four- 
velocity and Aa M is a four-vector containing NP CPT- 
violating coefficients. We approximate the 4-velocity of 
each B meson by the T(45) 4-velocity so that z is com- 
mon to each B in an event; in this way the 3- velocity 
direction becomes a function of the Earth position with 
respect to the Universe at the time of the T(45) decay, 
i.e. of the event's sidereal time, z is then expressed as: 

z ee z(£) = z +zi cos (Hi + (f>) (5) 

where Q = 2ir rad/sidereal day is the Earth rotation 
frequency (1 sidereal day ~ 0.997 solar days), t is the 
sidereal time, Zq and Zi contain NP coefficients. 

All the analyses presented in the following exploit the 
At dependence of the B-pair decay rates in order to ex- 
tract CP and CPT violation parameters, At being the 
difference between the decay times of the two mesons. 
The general form is: 
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In this equation Am (Ar) is the mass (lifetime) differ- 
ence between the B° physical states, while c+, c_ and 
s depend on the physical parameters \q/p\ and z, on the 
amplitudes of the particular processes under considera- 
tion, and the flavor of the B mesons at the time of decay. 



II. THE INCLUSIVE DILEPTON METHOD 
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Inclusive dilepton events, where both B mesons decay 
semileptonically (b — > X£v, with I = e or p,), represent 
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4% of all T(4S) — > BB decays and provide a very large 
sample with which to study T, CP and CPT violation. 
In the direct semileptonic neutral B decay, the flavor 
B° (B°) is tagged by the charge of the lepton £ + (£~). 
This study [6f is performed using a BABAR [7j data sam- 
ple equivalent to an integrated luminosity of 211 fb _1 
recorded at the T(AS) resonance (232 million BB pairs). 

Electrons and muons are selected using information 
from the tracking systems, the electromagnetic calorime- 
ter and the muon chambers. The track momentum is 
required to be between 0.8 and 2.3 GeV/c and certain 
dilepton mass ranges are vetoed in order to reject pho- 
ton conversions, J/?p and ip(2S) decays. The separation 
between direct leptons (b — > £) and background from the 
b — > c — ► I decay chain (cascade leptons) is achieved with 
a neural network that combines five discriminating vari- 
ables: the momenta and opening angle of the two lepton 
candidates, and the total visible energy and missing mo- 
mentum of the event. 

Since \q/p\ is expected to be small, we have de- 
termined the possible charge asymmetries induced by 
charge-dependent differences in the reconstruction and 
identification of electrons. The charge asymmetry of 
track reconstruction is measured in the data by compar- 
ing tracks reconstructed using only the silicon tracker 
with those passing the dilepton track selection, The lep- 
ton identification efficiencies are measured with a control 
sample of radiative Bhabha events for electrons. Based 
on tracking and identification, we fix the charge asym- 
metry for the direct electrons to — 1.2 x 1CP 3 . In 
the likelihood fit, we float the charge asymmetries 
and af, asc for direct and cascade muons. 

We model the contributions to our sample from BB 
decays using different categories of events each repre- 
sented by a probability density function (PDF) in At. 
The five categories are the following. First, the pure sig- 
nal events with two direct leptons are 81% of the BB 
events. Then, we consider two categories of cascade de- 
cays: those in which the direct lepton and the cascade 
lepton come from different B decays, and those in which 
the direct lepton and the cascade lepton stem from the 
same B decay. In addition, a small fraction of the dilep- 
ton events originate from the decay chain b — > t~ — > t~ 
which tags the B flavor correctly. Finally, the remaining 
events consist mainly of leptons from the decay of a char- 
monium resonance. The last component of the dilepton 
sample originates from mm-BB events, and has been esti- 
mated using off-resonance data. Figure[T]shows the result 
of the maximum likelihood fit in terms of the projected 
asymmetry Aqp/t between (£ + ,£ + ) and (£~,£~) dilep- 
tons and Aqpt between (£ + , £~ ) dileptons with At > 
and At < 0. 

For \q/p\, the most important systematic uncertainties 
are due to the correction of electron charge asymmetries, 
while the dominant systematic uncertainties on z are the 
imperfect knowledge of the absolute scale of the detector 
and the residual uncertainties in the tracker local align- 
ment. 
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FIG. 1: (a) A CP / T asymmetry between (£ + ,£ + ) and (£~ ,£~). 
(b) Acpt asymmetry between (£ + ,£~) dileptons with At > 
and At < 0. 



The results for CPT and CP violation parameters are 



\q/p\-l = (-0.8 ± 2.7 (stat .) ± 1.9 (syst . 
Imz = (-13.9±7.3( sta t.) ±3.2 (syst . 



) x 10~ 3 , 
x 10~ 3 , 



Ar-Rcz = (-7.1±3.9 (stat .)±2.0 (syst .)) x lO^ps- 1 . 

where Ar • Rez has been determined with the approx- 
imation Rez sinh(Ar At/2) ~ Rez • Ar • At/2 and z is 
considered as a constant parameter. 

A similar At fit is then performed, taking the sidereal- 
time dependence into account. We use the latitude and 
longitude of the BABAR detector, together with its ori- 
entation to determine the Lorentz boost direction of the 
T(4S) in the Earth frame. We convert the time record- 
ing of each event occurrence to calculate GMST (Green- 
wich Mean Sidereal Time) as specified by the U.S. Naval 
Observatory. Figure [5] shows the Aqpt asymmetry, inte- 
grated over the At variable, as a function of the GMST. 
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FIG. 2: The opposite-sign lepton asymmetry Acpt as a func- 
tion of GMST (in seconds) folded over a period of 24 sidereal 
hours. The curve is a projection from the two-dimensional 
likelihood fit onto the sidereal time axis. 



The results for sidereal-time dependent CPT violation 
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parameters are 
Imzo = (- 



-14.1±7.3 (Btat0 ±2.4 (ayirt0 ) x 10- 



Ar • Rcz = (-7.2 ± 4.1 (8tat0 ± 2.1 (syst . } ) x lO" 3 ps" 1 

Imzj = (-24.0±10.7 (Bta t.)±5.9 (syst .)) X 10" 
Ar ■ Rezi = (-18.8 ± 5.5 rstat .i ± 4.0^.0 x lO" 3 
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mesons which peak in the Ml distribution, like B + — > 
D**°£ + v e , D**° — > D*~tt + -. combinatorial events, corre- 
sponding to all BB decays not included in the previous 
categories; continuum events, coming from lepton-pair 
and light quark decays. Fractions of the different compo- 
nents in data are determined by means of fits to the Ml 
distribution: we obtain a signal fraction of (46.6 ± 0.9)% 
in the selected mass region — 4.0 < Ml < 2.0 GeV 2 /c 4 . 



III. THE PARTIAL RECONSTRUCTION 
METHOD 

In a complementary approach Q the partial recon- 
struction of the B° — > D*~£ + vi decay of one of the neu- 
tral B mesons is exploited. Though the available statis- 
tics is not as high as in the dilepton case, we can keep the 
charged B background at a lower level. At the same time, 
since the reconstructed and tag side are well defined, a 
procedure to determine particle detection asymmetries 
from data can be carried out. In this way, no dedicated 
control samples are needed to determine the asymmetry 
induced by the experimental cuts. A BABAR data sample 
equivalent to an integrated luminosity of 201 fb _1 at the 
T(AS) resonance (220 million BB pairs) is used for this 
study. 

Preliminary cuts for lepton selection and non-BB 
background rejection are similar to the inclusive dilep- 
ton analysis. Only the charged lepton from the B° decay 
and the soft pion (tt~) from the D*~ decay are recon- 
structed. The D° decay is not reconstructed, resulting in 
high selection efficiency. Due to the limited phase space 
available in the D* decay, we approximate the direction 
of the D* to be that of the ir s and estimate the energy 
Er>* of the D* as a linear function of the energy of the 
■k s - We define the square of the missing neutrino mass 
as: 



yfsjl-Eo, -FA -{PD*+P() 2 , (7) 
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where all quantities are defined in the T(AS) frame. We 
neglect the momentum of the B (approximately 0.34 
GeV/c), and identify the B° energy with half the total 
energy of the events (y/s/2). Eg and pi are the energy 
and momentum vector of the lepton and po* is the esti- 
mated momentum vector of the D*. The distribution of 
M^ is peaked for signal events, while it is spread over a 
wide range for background events. Backgrounds are also 
suppressed using cuts on the lepton momenta, on the soft 
pion momemtum, on the vertex fit probability, on |At| 
and a At- 

These selection criteria accept 470,877 events in the 
data sample (tagged events), and 5,291,868 partially re- 
constructed events that fail only the requirements of lep- 
ton tagging (untagged events). Fig. [3] shows the distri- 
butions of the squared neutrino mass for tagged events. 
Different components are recognized in the total sam- 
ple: signal events, due to D*~ ' l^vn events; peaking back- 
ground events, that correspond to decays of charged B 
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FIG. 3: Squared neutrino mass distributions for the tagged 
data events. The following fitted contributions are superim- 
posed: continuum (solid grey line), combinatorial (dotted), 
B ± peaking (dash-dotted), signal (dashed), all (solid black). 



As in the dilepton case, the contributions to our sample 
are modeled using probability density functions (PDF) 
in At, based on the signal-background categories defined 
above and on the true origin of the tagging lepton. 

An original method is used to determine possible 
charge detection asymmetries in the sample. We parame- 
terize the total number of expected events for each tagged 
and untagged category in a likelihood function, which 
contains a term coming from the shapes of the PDFs for 
tagged events, an extended term for the number of tagged 
events and an extended term for the number of untagged 
events. This method reduces the correlation between CP 
mixing asymmetry and particle detection asymmetry. 

Figure 2] shows the result of the maximum likelihood 
fit in terms of the projected asymmetries Acp/t between 
mixed positive and negative events, that we convert in a 
value for \q/p\. The main systematic uncertainties comes 
from the BB background and continuum charge asym- 
metries and the signal fraction determination from the 
M„ fit. The preliminary result is: 



\q/p\ - 1 = (6.5 ± 3.4 (stat , ± 2.0 (syst .)) x 10 



(8) 



which is compatible with the absence of CP violation in 
mixing at the 9.8% C.L. 



IV. CONCLUSIONS 

We have presented determinations of the parameters 
\q/p\, that is related to CP violation in B°B° mixing, 
and Imz, Ar ■ Rez, Imzi and Ar • Rezi, that measure 
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FIG. 4: Raw asymmetry between mixed positive (s r = 1, St = 
1) and mixed negative (s r = — 1, S* = —1) events with the 
PDF asymmetry derived from the fit superimposed. The 
asymmetry at low values of \At\ corresponds to the kaon- 
mistag contribution in decay-side tagged events. 

CPT violation in a theory that takes into account effects 
of the implied Lorentz violation. 

We exploit fits to At, the time difference between the 
two B decays in two complementary approaches. In the 
former, two high-momentum leptons are searched for in 
order to select inclusive semileptonic B° decays. In the 
latter, one of the B mesons is partially reconstructed in 
the semileptonic channel D*~t + v^, i.e. only the lepton 
and the soft pion from D*~ — » D tt~ decay are recon- 
structed, while the flavor of the other B is determined 
by means of lepton tagging. We use B°B° pair events, 
collected by the BABAR detector in the period 1999-2004. 

No evidence is found of CP violation in mixing with 
either of the two methods: both results are compatible 
with the SM expectations and with previously published 
BABAR results d,[l(|. The experimental values of A C p/r 
are used as constraints [llj , together with other measure- 
ments related to the unitarity triangle, in a global fit to 
the CKM theory parameters that also includes the pos- 
sibility of N P. Figure \5\ shows the most recent results of 
this fit in terms of the allowed regions for the parameters 
d and (j)d defined in Eq. (01, corresponding to a 95% 
confidence interval of [0.59,2.51] for C d and [-7.0,1.0]° 
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FIG. 5: Allowed regions for the NP parameters Cd and <f>d at 
the 68% (dark blue) and 95% (light blue) confidence levels. 

For the Imz, Ar • Rez pair, the results are compati- 
ble with the SM expectation (0,0) within 1.5o\ When 
sidereal-time dependence is taken into account, the re- 
sult in the Imzi-Ar-Reziplane has a significance of 2.2er, 
compatible with the absence of CPT violation in mixing 
(Figure [6]) . 




FIG. 6: Allowed regions for the CPT- violating parameters 
Im zi and Ar • Rezi at various confidence levels. The red star 
represents the SM expectation and the solid black ellipses 
correspond to 1, 2 and 3a significances. 
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